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Abstract 

A two-channel problem is considered within a method based on first order differential 
equations that are equivalent to the corresponding Schrodinger equation but are 
more convenient for dealing with resonant phenomena. Using these equations, it 
is possible to directly calculate the Jost matrix for practically any complex value 
of the energy. The spectral points (bound and resonant states) can therefore be 
located in a rigorous way, namely, as zeros of the Jost matrix determinant. When 
calculating the Jost matrix, the differential equations are solved and thus, at the 
same time, the wave function is obtained with the correct asymptotic behavior that 
is embedded in the solution analytically. The method offers very accurate way of 
calculating not only total widths of resonances but their partial widths as well. For 
each pole of the S'-matrix, its residue can be calculated rather accurately, which 
makes it possible to obtain the Mittag-Leffler type expansion of the S'-matrix as a 
sum of the singular terms (representing the resonances) and the background term 
(contour integral). As an example, the two-channel model by Noro and Taylor is 
considered. It is demonstrated how the contributions of individual resonance poles 
to the scattering cross section can be analyzed using the Mittag-Leffler expansion 
and the Argand plot technique. This example shows that even poles situated far 
away from the physical real axis may give significant contributions to the cross 
section. 
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1 Introduction 

Mankind has always been fascinated by striking phenomena such as thunderstorms, earth- 
quakes, and comets. As a part of mankind, scientists also are more attracted to drastically 
rather than smoothly varying quantities. Quantum physicists considering collisions of the 
molecular, atomic or sub-atomic particles, when find irregularities in their energy depen- 
dencies, try to assign certain quantum numbers to the outstanding features. These and 
other phenomena can often be regarded as resonances. Friedmann and Truhlar point out 
that "Probably the most striking phenomenon in the whole range of scattering experi- 
ments is the resonance" fl]. However, not all irregularities in the scattering data are true 
resonances [2]. 

Early studies of the energy dependence of nuclear scattering that led Bohr [3] to the 
idea of the compound nucleus (for more elaborate discussions see, for example, Ref. 
inspired similar ideas to appear in chemical physics. The notion of collision complex, used 
in atomic and molecular physics as well as chemistry, was developed in theoretical studies 
of the early 1970s ^i6j (see also Ref. [7]). 

However, it took about 20 years for indirect experimental evidence of the existence of 
such complexes to appear P E]. The first direct observation of an isolated resonance 
in chemical reactions only appeared around 2000 [71 |T0] . Still it is not absolutely clear 
that the results of Ref. fLO\ concerning the existence of the FHD collision complex in 
the F + HD — )■ FHD — )■ FH + D reaction, give an undoubted evidence. For example, 
in Ref.[TT] the doubts are expressed in terms of the definition of a resonance, which is 
rigorously defined as a quantum state that at large distances is a purely outgoing wave. 
Mathematically, this implies that the matrix element of the operator that transforms 
the incoming wave into the outgoing wave (the S'-operator) , has a pole at the resonance 
(complex valued) energy [12], [13] . 

A theoretical analysis of intermediate states, in the collisions like F + HD, will eventually 
need more than one electronic potential energy surface, a feature which is not easily 
achieved [14] . It is thus, from a theoretical point of view, reasonable to begin with atom 
ion collisions, which can be represented by a set of coupled potential energy curves. 
Looking back in the literature, we find that a few resonances were found in a theoretical 
study of the N^~^ + H ^ iV^"*" -|- reaction by McCarroU and Valrion[15j. Barany et 
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al. [15] and Rittby et al. [H] extended this work and could assign complex eigenenergies 
and their rovibrational quantum numbers to not only the resonances found in [TB] . These 
resonances were, as in the work by Elander and coworkers [T8t [T9 | [20l [2T] . forming a sting 
that started with a narrow resonance and propagated out in the complex energy plane. 
The group lead by Shimakura, studied several charge-transfer reactions with the initial 
channel N^+ + H (see Refs. [221 [221 121125]) while Bacchys-Montabonel and Ceyzeriat [25] 
did calculations for the ion-atom collision Si'^'^ + He. Both groups predicted resonances 
in these reactions. Experimentalists at Oak Ridge made an effort to find evidence for 
these type of resonances [271 12H] but the obtainable energy range and resolution did not 
permit their detection. The conclusion was that due to extremely high Coulomb barrier 
the relative energy (~meV) at which these resonances appear, could never be achieved in 
collision experiments. 

The interest to the resonant charge-transfers and other similar processes did not wane. 
The new idea is to observe the collisions between particles moving in the same beam, so 
that their relative kinetic energy is very low. Electrostatic storage rings for this purpose 
are currently being built. For example, the Desiree project in Stockholm [29] is planned 
to have the possibility of merging a stored multiply charged ion beam with a beam of 
neutral atoms. Similar experiments are also planned at other facilities. 
In the theoretical studies, the most ubiquitous formula associated with resonance phe- 
nomena, is the one bearing the names of Breit and Wigner. It was invented for phe- 
nomenological description of the peaks in the neutron scattering cross section [30j and is 
valid in the neighborhood of an isolated resonance. From textbooks, we learn that if a 
resonance is embedded in a direct scattering continuum, the cross section profile does not 
have a Lorentian shape as predicted by the Breit- Wigner formula, but should rather be 
described by the more general Butler- Fano parameters [3T] . 

The original ideas of Breit and Wigner [30] and Fano [31] may make us believe that reso- 
nances only give a signature to a cross section as peaks or similar irregularities. However, 
model potential studies of the Schrodinger spectrum [221 [22] found that the S-matrix 
poles associated with resonances, form long (most likely infinite) strings in the complex 
energy plane. The influence of these poles on the cross section could be more complicated 
than simple generating of the peaks or bumps. The resonances may interfere with each 
other and thus cause an intricate energy dependence of the cross section. 
The most controversial example of such intricate interference of resonance states is the 
problem of barrier states in chemical physics. According to Friedmann and Truhlar, 
"metastable states associated with barriers, are associated with poles of the S-matrix just 
as definitely as are trapped states associated with standing waves in a well" and these 
metastable states are the bottle necks for chemical reactions (see Refs. [1] and [31]). 
In relation to barrier scattering reactions like H + HD — ?■ D + H2, these ideas are dis- 
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cussed by Althrope et al. where they report "an intriguing forward-scattering peak" 
observed at a coUision energy ~1.64eV in a state-to-state differential cross section. From 
the point of view of Friedmann and Truhlar [Tj, this observation could be a manifesta- 
tion of quantum barrier states suggested by them. This is opposed by Manolopoulos |36] 
who claims that isolated (i.e. true) resonances and barrier type resonances are "mathe- 
matically different and the physical implications of this mathematical difference are too 
significant to regard the two situations as a manifestation of the same phenomenon" [55] . 
Other similar results on "quantized bottleneck state(s)" are reported in Ref. [37]. It 
should also be mentioned that the quantized bottleneck ideas do have implications on the 
understanding of the so called transition states discussed in chemical physics ( See Ref. 
[55] for a review). 

The controversies, like the barrier states, can only be resolved if we reach a clear un- 
derstanding of how much a given isolated or barrier type resonance influences the cross 
section. The importance of this is emphasized by the authors of Refs. [H |3ll [351 ISHl |37] 
and many others. A number of methods were developed for the purpose of analyzing the 
roles of individual resonances in the scattering cross section. A comparison between the 
Mittag-Leffler expansion approach (see, for example, Ref.[TB]), the Siegert pseudo-state 
method [39J, and the complex Kohn variational method was recently published in Ref. 
[10] . Based on the Titchmarsh- Weyl theory, the authors of Ref. [31] developed a method 
for decomposing the scattering information into the resonance and background contribu- 
tions. One of us (N.E.) has earlier analyzed the Mittag-Leffler method to understand its 
computational properties [IHl IMl [21] • In Ref. [21], we also investigated the influence of a 
string of resonance poles yielding somewhat surprising results. 

One drawback of the Mittag-Leffler expansion method is that it may be hard to apply 
it to the systems involving more than two particles in at least one of the channels. In 
this respect the complex variational Kohn method, as described by Nuttal and Cohen 
|42] . Rescigno and McCurdy[13], and recently in the review by Moiseyev [H], has an 
advantage. It can be formulated in terms of matrix elements of the interaction potential 
and the complex eigenvalues of the corresponding analytically continued Hamiltonian 
[44] . In this way the Kohn method is applicable whenever a basis set representation of 
the problem is available. It can therefore be extended to the 3-D problems such as, for 
example, the colhsions complex FHD etc. [HI [SI HD]- However, being based on the basis 
expansion, the complex variational Kohn method has its disadvantages as well. These 
are related to the problem of convergence and completeness of the expansion. It is thus 
of importance to compare and analyze the three above mentioned methods in parallel to 
understand and utilize their respective powers in analyzing and predicting experimental 
results. 

As a prerequisite for the analysis of the cross sections of chemical reactions, one needs 
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a reliable method for solving the multi-channel Schrodinger problem. As we already 
mentioned, in theoretical analysis of atomic complexes, like F + HD, it is very common 
to represent the atom-ion interaction by a set of coupled potential energy curves. 
In this paper, we adopt a mathematically rigorous definition of a resonance, namely, 
as a point on the unphysical sheet of complex energy surface at which the Jost-matrix 
determinant is zero. At such points, the S'-matrix has poles since this determinant is 
in the denominators of all its matrix elements. The so-called redundant (non-resonance) 
poles that the S'-matrix may also have (see, for example, Ref.[l5]), are thus avoided 
since they are associated with singularities of the numerators of its matrix elements. 
This mathematically-inclined definition is consistent with physics in that at thus defined 
resonance points the wave function has the Gamow-Siegert asymptotics (pure outgoing 
waves), but it does not say anything about observable irregularities of the cross section. 
Therefore the term "resonance" has here a more general meaning not limited to long-lived 
(narrow) states. 

In the paper, we describe a mathematically rigorous and numerically accurate and stable 
method for solving the two-channel (generally, A^-channel) quantum mechanical problem. 
For the last ten years various aspects of this method were developed in Refs. [IHl SZl 
SSlllSlEniEIllSiESlIMllSSlESlEZlEB]. Here we suggest an approach to analyzing 
the contribution of resonance S'-matrix poles to the scattering picture. Our approach 
is based on a Mittag-Lefiler type expansion of the S'-matrix, for which we calculate the 
residues of the S'-matrix at the resonance poles and its Cauchy-type contour integral. 
The roles played by individual resonances are deduced using the Argand plot technique. 
Practical application of the method is demonstrated by the example of the Noro and 
Taylor two-channel model |59] . 

2 Two-channel Schrodinger equation 

Let us consider a quantum mechanical two-body problem which, after separation of the 
motion of its center of mass, is reduced to an effective problem of one body whose dynamics 
is governed by the Hamiltonian 

H = Ho + U + h, (1) 

where the terms on the right hand side describe its free motion {Hq), the interaction forces 
(f/), and the internal dynamics in the body [h), respectively. Since our main objective is 
to study how the resonance poles of the S-matrix manifest themselves in experimentally 
observable cross sections, we avoid unnecessary complications associated with spins of the 
particles and with long-range Coulomb forces. Furthermore, for the sake of clarity, we 
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consider here a two-channel system, although all the equations can be easily generalized 
for an arbitrary number of channels (see Ref. as well as for charged particles with 

non-zero spins (see Refs. 150]). 

Therefore we assume that the colliding particles are neutral and spinless, and the Hamil- 
tonian h acts in the space spanned by only two of its eigenstates, 



h\n) = En\n) 



n 



1,2 



(2) 



1. e. 



/i = Ei|l)(l| + E2|2)(2| . (3) 

In its turn, this implies that the operators Hq and U are 2x2 matrices in the subspace 
associated with the internal degrees of freedom. Therefore the total Hamiltonian ([T]) taken 
in the coordinate representation and sandwiched between {n\ and has the following 
matrix representation 



/ ^2 



H 



2/i: 



■A,-+f/n(f) + Ei 



f/i2(rl 



t/2i(r-0 



2/i; 



■A,^+f/22(r) + ^2 



(4) 



where r is the relative coordinate and the subscripts label the channels. A solution of the 
corresponding Schrodinger equation 



Hi) = Ei) 



(5) 



is a column matrix 



composed of the f-dependent coefficients in the expansion of the state vector over the 
channel (internal) states, i.e. 



(6) 



Continuing to keep unnecessary complications away, we assume that the interaction po- 
tentials are spherically symmetric which means that the orbital angular momentum £ 
associated with r, is conserving and hence is the same for the in and out channels. Fur- 
thermore, for the matrix elements of the potential, we assume that 



\Unr,'(r)\rdr < oo 



(7) 
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i.e. that their possible singularity at r = is less than and they vanish at infinity 
faster than Expanding the channel wave function in partial waves, 



ipn{E,r) = Yem{r) 



(8) 



we obtain a set of coupled radial Schrodinger equations 

nu = Vu , (9) 
where m is a column matrix consisting of Ui{E, r) and U2{E, r), matrix elements of V are 



K.n'lr) 



■U„n'(r) 



and 



n 



\ 



V 







with the channel momenta kr, defined as 



{E-E„ 



+ kl- 



+ 1) 



(10) 



(12) 



As is well known (from the general theory of ordinary differential equations), a coupled 
set of N equations of the type ^ has 2N linearly independent column solutions, and 
only half of them are regular at r = 0. Since in our case = 2, we have two independent 
regular column-solutions that can be combined in a square 2x2 matrix 



0(E,r) 




k2{E,r) 
h2{E,r) 



(13) 



where the first subscript labels the channels and the second is the solution number. For 
the independent regular solutions, we use the symbol to distinguish them from the 
physical solutions u. By definition, they obey the same equation, namely. 



'H(p = V^ . 



(14) 



Since all physical solutions u{E,r) must be regular at r = 0, anyone of them is a linear 
combination of the regular columns, i.e. 



(15) 





or in a more compact form 

/ Ml \ / 011 012 \ / Cl \ 

(16) 




^21 




The combination coefficients Ci and C2 should be chosen in such a way that guarantees a 
given (physical) behavior of the solution when r — )■ oo. 



3 Transformation to first-order equations 

If the right hand side of Eq. flTl]) were zero, the equations of the set would be uncoupled 
(independent) and could be solved analytically. Indeed, the Riccati-Hankel functions 
h^f^\knr) are solutions of the equation 



ht\Kr)=0 



Therefore the matrices 



, , , h'f \k,r) 

wt\E,r)={ I (17) 

h'f\k2r) 



and 



h['^\kir) 

h[~^\k2r) 
both solve Eq. (fT4|) when its right hand side is zero, i.e. 

W/"/°"*^(E,r) = . (19) 

These matrices represent the incoming and outgoing spherical waves which form the 
asymptotics of the solution at large distances when the potential vanishes. In order 
to guarantee correct asymptotic behavior of the solutions (fT3!) of Eq. (IT^ . let us look for 
them in the following form 

<j){E, r) = Wt^ {E, r) {E, r) + I^f {E, r)^(°"*) {E, r) , (20) 

where J^("V°"*)(E, r) are new unknown matrix functions. In the theory of ordinary dif- 
ferential equations, this way of finding solution is known as the variation parameters 
method. 
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Since we replaced one unknown matrix with two unknown matrices ^ they cannot 

be independent. We therefore can impose one arbitrary condition that relates them to 
each other. As such condition, we choose the following equation 



(Jj I (Jj I 



(21) 



which is standard in the variation parameters method and is called the Lagrange condition. 
It implies that, when calculating first derivative of (j){E,r), we only need to differentiate 
the spherical waves. 



d 



j-(t>{E,r) = W',^''\E,r)J^^-^\E,r) + r) J-^""*^^, 

(Jj f 



(22) 



Let us substitute the ansatz (l20i) into Eq. ( |T^ and thus obtain the corresponding equa- 
tions that determine J-'^™) and The second derivative of (j){E^r) is 

j2 

^ p^"(in)jr(m) _^ y^"(out)^(out) _^ yj/'Mj-'(in) _^ y^' (out) (out) _ 

This means that 

where the first two terms disappear in accordance with Eq. ( |T9l) . Therefore Eq. 
takes the form 



,(in)^,(i,)^ 



/(out) J-/ (out) 



V 



(in)^(in) ^ w^-Vjr 



(out) -77(out) 



From the Lagrange condition (ETj), it follows that 



-pi (out) 



W, 



(out) 



-1 



W(i-)jr'(in) _ 



(23) 



(24) 



It should be noted that since matrices VF^^™^""*-* and their derivatives all are diagonal. 



they commute with each other. Substituting Eq. ( I24l) into Eq. ( I23l) and using this 
commutative property, we obtain 



p (in) 



^(out)^/(in)_^/(out) (in) 



y^(in)jp(in) ^ p^(out)^(out) 



Similarly, from the Lagrange condition ( 1211) . it follows that 



jrl (in) 



(in) 



-1 



y^(out)jr' (out) 



(25) 



(26) 
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and therefore Eq. 

-pi (out) 



gives 



^(in)^,(out) 



^/(in)^(out) 



-1 



(out) -77(out) 



(27) 



The first factors in Eqs. fl2^ and fl27|) differ only by the sign and actually are equal to the 
inverted Wronskian of the incoming and outgoing spherical waves, which can be found in 
explicit form. Indeed, knowing the Wronskian of the Riccati-Hankel functions, 

h^[\kr)^h^^\kr) - ^h^[\kr)hf\kr) = 2ik , (28) 



we obtain 
where 



^an)^^'(out) 



dr 



^,(in)^Kt) 



2iK 



K 



(29) 
(30) 



h 

k2 

Therefore the Schrodinger equation is equivalent to the following set of coupled 
equations of the first order 



dr 

_^jr(out) 
dr 



1 



2i ^ 



2i ^ 



(in)^(in) ^ wvou.;^, 



(out) ■77(out) 



y^(in)jp(in) ^ y^(out)^(out) 



(31) 



Boundary conditions at r = for equations flHTl) follow from the requirement that the 
solution fl2U]) must be regular. This can only be achieved if J-'^^^\E,r) and J-'^°^^\E,r) 
become identical when r — )■ because in such a case the singular parts of h^f \kr) and 
h^^\kr) cancel each other, i.e. 



hY\kr) + h'f^>{kr) = 2je{kr) 



Here the Riccati-Bessel function ji is regular, namely. 



je{kr) 



TT 



-{kr) 



(32) 



(33) 



r^o 2^+ir(£ + 3/2) 

Matrices J^^™\E^r) and J^^°^^\E^r) must become diagonal near r = in order to guar- 
antee linear independence of the columns of (j){E,r). Thus, the following boundary con- 



ditions 



hm:F^:/-^{E,r)=6„ 

r— >0 



(34) 



are appropriate since we have the freedom to choose the normalization of the solution. 
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4 Jost matrices 



The differential equations (13T|) can be numerically solved from the origin to a sufficiently 
far point r = R where the potential vanishes (causing the right-hand sides of the equations 
to disappear) and therefore J-'^^^^°^^\E,r) become constant. These constants 

= lim J^('°)(E,r) (35) 

r—^oo 

and 

r—^oo 

are the Jost matrices that determine asymptotic behavior of the fundamental system of 
regular solutions 

(j){E,r) wt\E,r)f^'''\E) + wt"'\E,r)f^°'''\E) . (37) 

It is worthwhile to mention that we use here the notation /('°/°"*)(_E) which is different 
from the traditional notation such as f^'^\±k). There are two reasons for this. First of all, 
we do not fix the normalization of the regular solution 4>{E, r). As a result, both f^^^\E) 
and /(°"*)(E) can have an arbitrary common factor. We are not concerned with this 
factor because no observable quantity depends on it. In contrast to a majority of other 
studies, we leave the normalization of the Jost matrices free, and therefore need a notation 
that is different from the traditional. The second reason is that the superscripts {in) and 
(out) are unambiguous and thus we avoid possible confusion caused by the existence of 
notations with opposite signs for the same Jost matrices. 

The regular solution 0(-E', r) consists of two linearly independent columns which constitute 
the fundamental system of solutions. Any other solution is a linear combination of these 
independent columns. The physical wave function fllSp is one of such combinations, 

u{E,r) = (j){E,r)c , (38) 

where 

u = ( I and c = I ] . 

\U2 J \C2 J 

At large distances, we therefore have 

u{E,r) wt\E,r)f^''\E) f + iy/°"*)(E, r)/(°"*)(E) f • (39) 

■r-s>oo V '^2 / V '^2 / 
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4.1 Spectral points 



There are certain discrete points in the complex ii^-plane (bound, resonant, and virtual 
states), at which the physical wave function has only outgoing waves in its asymptotic 
behavior. At these so called spectral points the combination coefficients ci and C2 are such 
that 

This homogeneous equation has a non-trivial solution if and only if 

det/('°)(^) = . (40) 

As is seen, the distribution of the spectral points in the complex i?-plane (determined 
by this equation) does not, as was mentioned earlier, depend on the normalization of the 
Jost matrix. 




4.2 Scattering 

For a scattering state, the wave function has both incoming and outgoing waves at large 
r. If An and are the amplitudes of the incoming and outgoing (scattered) waves in 
the n-th channel, then the scattering boundary condition (at r — )■ oo) reads 

h[-\hr)A, + h[-^\hr)B, 
u{E,r) I I (41) 

h^f\k2r)A2 + h^^\hr)B2 

= wt\E,r) 1^ j +wt''\E,r) I • (42) 

Comparing Eqs. and (15^ . we see that 

^1 \ i-(in)/ r7'\ / Cl 



and 



J - K; 1 . (44) 
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This implies that the amphtudes of the incident and scattered waves are related by a 
2 X 2-matrix , S{E) as 

; S ) - it)- 

This S'-matrix 

5(E) = (46) 

does not depend on the choice of An. We can therefore clarify the physical meaning of 
its matrix elements by considering special cases with simple choices of Ai and A2. If 
the incident wave is purely channel 1, i.e. Ai ^ and A2 = 0, then Rn = Bi/Ai and 
T12 = -82/^1 are the amplitudes of the probability that the incoming wave is returning 
in the channel 1 and is transmitted into channel 2, respectively. Similarly, ii Ai = 
and ^2 7^ then R22 = -B2M2 and T21 = B1/A2 are the amplitudes of the probability 
that the incoming wave is returning in the channel 2 and is transmitted into channel 1, 
respectively. Substituting = or A2 = into the equation 



Bl\( SuA^ + ^12^2 
B2 / \ 5*21^1 + 5*22^2 



(47) 



we see that the S'-matrix consists of these transmission and reflection amplitudes, namely, 

^^^^ - [ T,2{E) , R22{E) ) ■ ^48) 

In general the energy E is complex and only for real energies the quantities Rnn and Tnn' 
have the simple physical meaning. 

Also note that the S-matrix and therefore the observable transmission and reflection 
amplitudes do not depend on the normalization of the Jost matrices. Indeed, any common 
factor of Z*-™'' and /(°"*) cancels out in Eq. ( H6|l which supports our suggestion to move 
away from their fixed normalization. 



4.3 Partial widths 

At the energy points corresponding to bound states, and at complex energies in the 
resonance region (fourth quadrant) of the i?-plane the elements of the S'-matrix have 
special properties as well. From Eqs. f HBj) and (HUj) . we see that at every bound state and 
resonance 

E,,, = E,-i^ (49) 
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the S'-matrix has a pole. If the energy is close to Ej-^s then 



Snn'iE) ^ const (l- ^ 'VrnTn' \ ^ ^ 

where Fi and r2 are the partial widths such that together they form the total width 



r = Ti + r2 



(51) 



of the resonance, and the ratios Fi/F and r2/r are the probabilities that the resonance 
will decay through (or can be excited from) the first or second channels, respectively. As 
it follows from Eq. ( I50l) . the partial widths r„ can be found as the limits 

r„ = lim {E - Er^s)Snn{E) . 

However, in numerical calculations, finding a limit of a singular functions is not an easy 
task. A way to avoid this difficulty was suggested by Masui et al. in Ref. [SI]- From Eq. 
fISU]) it is clear that 



lim 

E^E res 



SuiE) 



£1 

r2 



(52) 



S22iE) 

Therefore the partial widths can be found with the help of Eqs. (!5T!) and ( |52|) . if F is 
known and if there is a procedure for the S'-matrix calculation at complex energies. 
We can significantly simplify the task of finding the limit ( l52l) by using Eq. (H6i) . Indeed, 
both 5*11 and 6*22 have singularities at E^es because of vanishing determinant of Z^™-*, which 
cancels out in their ratio. Making explicit inversion of the 2x2 matrix Z*^™) and using 
Eq. f H6|) . we obtain 



£1 

r2 



/•(out) r{in) /.^^ouL; r 

In J22 " J12 J'. 



(out) ('(in) 



21 



f (out) f (in) 
i22 ill 



f'. 



(out) (-(in) 
/l2 



21 



(53) 



E — Eri 



where no singularities are present. 



5 Complex energies 

The differential equations fl3T|) enable us to obtain a complete solution of the two-channel 
(generally, iV-channel) problem at any energy of physical interest. Their advantage over 
the corresponding Schrodinger equation becomes especially evident when we consider 
complex values of the energy. 
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5.1 Riemann surface 



Since the functions J-'(in/out)(£', r), and therefore the Jost matrices , depend on the en- 
ergy E via the channel momenta ki = a/ {2jji/h?){E — Ei) and k2 = a/ {2jj2/K^){E — E2), 
there are two square-root branching points for these matrices in the energy plane : E = Ei 
and E = E2. This means that if we make two full circles around either of these points 
then matrix elements of return to the same values from which the circling was 

started (one full circle is not enough). This is because the channel momentum 

kn = ^/{2fiJh^)\E-En\ex^ = V(2/i„//l2)|E-E„|e*^"/2 (54) 

comes to its initial value when Xn — > Xn + 47r, where \E — En\ and Xn are the polar 
coordinates of the point E on the energy-plane relative to the branching point E^, i.e. 

E = En + \E - Er^le'^" . (55) 

In other words, elements of the Jost matrix have two different values at each point E on 
this circle. 

In order to make the Jost matrix a single-valued function of E, we can assume (as is usual 
in the complex analysis) that the complex energy forms the so-called Riemann surface 
consisting of several parallel sheets. When doing the first circle around a branching point, 
we are moving on the first sheet and then continue on the second one until coming back 
to the first sheet after completing the full two circles. Such continuous transition from 
one sheet to another is possible if we make a cut from the branching point to infinity, and 
connect opposite rims of the cuts on the two sheets (see Fig. [1]). 

As is usual in the scattering theory, we make straight line cuts from both branching 
points to infinity along the positive real axis. Each of the two sheets related by the 
first branching point, is further branched at the second branching point. Therefore the 
full Riemann surface consists of four parallel sheets. We can reach any of these sheets 
by making an appropriate number of circles around the first and the second branching 
points. 

The physical energy (at which the scattering takes place) is on the positive real axis. 
We choose the cuts and their interconnections in such a way that the physical scattering 
energies lie on the upper rims of the both cuts. Starting from these physical energies 
and moving in the anti-clockwise direction around all branching points, we cover the so- 
called physical sheet of the energy plane. According to Eq. f l5^ the channel momenta 
corresponding to this sheet, have positive imaginary parts. The bound states are also on 
the physical sheet (this is necessary to guarantee the exponential attenuation of their wave 
functions). The resonances, however, correspond to zeros of the Jost-matrix determinant 
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at the momenta with negative imaginary parts and therefore he on the unphysical sheet of 
the Riemann surface. The resonances that are able to decay into both channels must have 
Im (ki) < and Im (^2) < 0. This means that they lie on the sheet which is unphysical 
with respect to the both branching points. 

In practical calculations, we can choose point on one of the four Riemann sheets by 
selecting appropriate signs in front of the square roots for ki and k2, i.e. appropriate 
signs of their imaginary parts. These signs must be such that the imaginary part of the 
momentum is positive for a closed channel and negative for an open one. This corresponds 
to the choice between the physical and unphysical sheets of the Riemann i?-surface and 
ensures correct asymptotic behavior of the wave function. 

5.2 Complex rotation 

Similarly to the single-channel case (see Refs. SH]), the scheme described in previous 
sections, can be easily implemented only for bound and scattering states. In the resonance 
domain of the complex ii^-plane (below the real axis) the Riccati-Hankel function h^^\kr), 
as can be seen from its asymptotics 



diverges when r — t- 00. As a result the right hand side of the first matrix equation of the 
set fl3T|) diverges and hence the limit fl35|) cannot be calculated. It should be emphasized 
that this fact does not mean that this limit does not exist. It do exist, but simply moving 
along the real r-axis, we cannot reach it. It is easy to see, that 



When Im kr = this function remains finite (oscillates) at large r. The condition Im fc^r > 
for asymptotic vanishing of the function h^^\T]n, knr) involves the imaginary part of the 
product knT but not of the momentum alone. This offers an elegant way to extend the 
domain of the E'-plane where the limit (135|) can be calculated, to practically whole E- 
plane. Indeed, if, for example, Im /c„r is negative we can always make it positive by using 
complex values of r. This of course requires that the potential is defined for complex r 
and tends to zero when |r| — )■ 00 at least in certain sector of the complex r-plane. We 
assume that it vanishes faster than when r — )■ 00 along any line 




(56) 




|A;r|-s>0 



(57) 



r = zexp[v 



z>0, 



(58) 
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for the rotation angle 9 in the interval < |^^| < ^^max < ^/2. Then, firstly, according to the 
existence theorem, the solutions 0„„'(i?,r) of the Schrodinger equation are holomorphic 
functions of r within the corresponding domain of the r-plane and, secondly, they have 
the asymptotic behavior of the type fl37|) along any line (|58|1 . Moreover, the coefficients 
j(m/out)^^^ in such asymptotics are the same for all choices of the rotation angle because 
they do not depend on r. 
Similarly to Eq. f l57p . we have 



Therefore, the limiting values oo) can be found as the corresponding coeffi- 

cients in the asymptotics of the functions ^e*^). It can be shown that with 

Im fcr > the right hand side of the first equation of the set (131!) at large distances tends 
to zero, which means vanishing of the derivative on the left hand side, i.e. that J-"*^'") 
becomes constant (reaches its limit). Similarly, it can be shown that the limit fl36p can be 
reached if Imfcr < 0. Both limits, i.e. f^\E) and can be found simultaneously 

only when Im fcr = and at all spectral points. The line Im fcr = therefore serves as the 
dividing line that separates two domains of the complex energy plane. The Jost matrix 
f'y™)i^E) can be calculated above this line, while f^°"^\E) below it. 

By considering complex r, we actually do the analytic continuation of f'^™/'^^^)i^E) across 
the dividing lines (the unitary cuts) to the domains where Eqs. ( 13T|) do not give finite 
values for these matrices. The dividing lines can be turned down wards to expose the 
resonance spectral points, by rotating r as given by Eq. (l58!l . Indeed, if Xn is the polar 
angle parametrizing the position of a point on the i?-plane relative to the branching point 
En, then by choosing large enough 9 we can make ImA;„r, 



positive even when Xn is negative (when the point E is below the cut, i.e. on the unphysical 
sheet) and vice versa. From the last equation is clear that when 6* > both dividing lines 
are turned down by 29 (see Fig. [2]). 

The fact that f^^^\E) and f^°^^^(E) can be found simultaneously only on the lines corre- 
sponding to the unitary cuts and at the spectral points, does not pose a problem. Indeed, 
for bound states Imfcr > and we can locate them using Eq. fHU]) without the rota- 
tion. As soon as a bound state has been located at a spectral point Ef,, we can calculate 
both J-'^'^\Eb,r) at this point and thus obtain the wave function. The scattering takes 
place exactly on the real axis that coincides with the unitary cut where both f^^^\E) 
and /(°"*)(E) can be calculated without rotation. To locate resonances, we turn the cuts 




(59) 



ImknT = Im (|A;„|ze^(^+'^"/2)) = \ki\zsm{9 + x„/2) , 



(60) 
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down opening the unphysical energy sheets, and again use Eq. fl4UI) which involves only 
f(^^)(^E). If we need to obtain the partial widths (which requires the knowledge of the 
f('^^'^)(^E) as well), we have to repeat the integration of the differential equations at the 
resonance energy with such rotation angle that the resonance in question is below the cut. 
It should be also mentioned that the border separating the two domains of the complex 
ii^-plane is a line only in the case of potentials slowly vanishing at infinity (for example, 
as 1/r"). If, however, the potential decays exponentially, then f^™\E) can be found not 
only above the dividing line but also within a band below this line and f^"^'^\E) within 
a symmetrical band above the line (see Ref. [50]). The faster the potential decays the 
wider this band is. 



5.3 Integration path 

As was shown above, when solving the differential equations fl3T|) . we have to move from 
the origin to infinity (actually to a large |r| = R) along the ray (|58|) with generally 
nonzero 6. In numerical calculations, however, this straightforward approach is not always 
practical. As was mentioned earlier, in the immediate vicinity of the point r = 0, the 
singularities of the Riccati-Hankel functions cancel each other in accordance with Eq. 
( 132|) . This may make Eqs. ( 13T|) numerically unstable near r = 0. Such an instability is 
further aggravated when the angular momentum i becomes large. 

This problem can be easily circumvented in two ways. One way is to solve the matrix 
Schrodinger equation f lT^ on an interval from a sufficiently small r = r^^^ to an inter- 
mediate point r = b. Having thus calculated the matrices (j){E,b) and (j)'{E,b), we can 
obtain the values J'^*"/""*) (E, b) using the relations 

^(in/out) = ±Yi^'' [iy;^°"'/'"^(E,r)0(E,r) - iyj°"*/'"^(E,r)0'(E,r) 



which follow from Eqs. (1201) and ( 1221) . Then we can turn to the ray (15811 along the path 
shown in Fig. |3l and safely proceed with Eqs. (|3T|) to a sufficiently large z = R, which 
will serve as the "infinity" . 

Another way is to rearrange Eqs. (|3T1) on the interval [rmim b] in order to avoid the 
singularity cancellations. Since {h^(^^ + h\ ^)/2 = j£ and (/i^^"* — h\ ^)/2i = ne, where is 
the Riccati- Neumann function (which is singular at r = 0), we may introduce a new pair 
of matrices 

AiE, r) = J^("^) {E, r) + {E, r) (61) 

and 

B{E, r) = % [J^('") [E, r) - J^(°"*) (E, r)] , (62) 
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which transform the ansatz (I2UI) into the form 



r) = ME, r)A{E, r) - J\f,{E, r)B{E, r) , (63) 

where 

/ j,{hr) \ 
ME,r)=[ (64) 

and 

/ neihr) \ 
J^i{E,r)= . (65) 

V ni{k2r) J 

The corresponding hnear combination of Eqs. 0311) gives their alternative form 

d 



-A = -K-^UtV\JiA-UiB] 
dr 

= -K-'jfV[J,A-M,B] 

dr 



(66) 



with the boundary conditions that follow from ([3 



lim^„„/(£;,r) = 25nn' , lira Bnn'{E,r) = . (67) 

r— >-0 9 — ^0 

The representation of (j){E,r) in terms of A, B and j^(™/out) jg equivalent. However, from 
a practical point of view, near the origin it is more convenient to use Eqs. (166|1 which 
do not have the problem of singularity cancellation and thus are numerically stable. At 
large distances, we have to solve Eqs. in order to obtain the Jost matrices {A and 
B involving both generally do not converge with any rotation angle). Thus, 

we can follow the same path shown in Fig. |3l making transformation from {A, B} to 
|jr{in)^ jr(out)j at r = 6. 

In principle, from the very beginning when solving either the Schrodinger equation or 
Eqs. ( 166|) . we could move along the ray ( l58l) . This however could be a source of numerical 
instability of another kind. If the potential has exponential functions, it is oscillating with 
complex r. Moving along the real axis allows us to avoid this problem. Moreover, the 
potential may be given in an analytical form only for large r while in the inner region it 
is known in the form of a table along the real axis. 
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6 Numerical example 



In order to demonstrate the efficiency and accuracy of the proposed method, we apply it 
here to a particular two-channel problem. Another purpose of the calculations reported 
further down, is to learn how individual resonances contribute to observable scattering 
cross sections. 

6.1 Potential 

As a testing ground, we chose the model proposed by Noro And Taylor [59] and used in 
many other publications since then. Their two-channel potential 

/ -1.0 -7.5 \ 

U{r) = r^e-' , (68) 

V -7.5 7.5 / 

is of a short range and apparently obeys the condition ([7]). The units in this model 
are chosen in such a way that /ii = = = 1. In the tables and figures given in 
the subsequent Sections, the units for the energies and cross sections are therefore not 
indicated. 

The threshold energies are Ei = and E2 = 0.1. The potential fl68|) has an attractive 
well in the lower channel, a barrier in the upper channel, and rather strong coupling 
between the two channels. This means that it generates a rich spectrum and represents 
a real challenge to any method that is designed to analyze the roles played by individual 
resonances in the whole picture of the scattering process. 

All the calculations presented here, were done for i = 0. This is not because our method 
has some limitations or needs simplifications. It can handle any value of the angular 
momentum (even complex). Since however this is an artificial model and there are no 
partial waves that would present a special physical interest, we chose the lowest one. 
Moreover, to the best of our knowledge, in all other publications only the S'-wave scattering 
was considered and we therefore have to do the same to be able to compare some results. 

6.2 Spectral points and cross sections 

Using the Runge-Kutta-Fehlberg method [60], we solved Eqs. ( 166|) from rmm = 0.0001 
to 6 = 1 along the real axis. Then with the same numerical method, we solved Eqs. 
(pTj) from r = 6 to |r| = 30 along the path shown in Fig. |3l For searching the bound 
states and calculating the S'-matrix for real energies, the rotation angle 6 was taken to be 
zero. When the resonances were located and the S'-matrix was calculated in the complex 
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i?-plane, we used appropriate values for this angle (up to OAtt) that turned the unitary 
cuts (see. Fig. [2]) far enough to open necessary domains of the complex plane. 
The spectral points nearest to the origin of the E-plane are given in Table [1] where 
the corresponding results from some other publications are presented for the sake of 
comparison. The S-matrix poles corresponding to these spectral points are shown in 
Fig. HI The potential fIBS]) supports four bound states and generates a string of nine 
resonances at energies above the second threshold. All other spectral points are sub- 
threshold resonances and lie at negative energies. 

The two lowest bound states lie below the minimum of the attractive potential in the first 
channel. This shows that the system under consideration is essentially diabatic as one 
should expect with strong cross-channel coupling in the potential fl68p . 
The S'-wave scattering cross sections for the transitions from channel n to channel n', 



presented in Figs. [5l El and [71 show several bumps and peaks. Our traditional under- 
standing of quantum scattering, preoccupied with the Breit-Wigner picture, would push 
us towards associating these cross section irregularities with individual resonance poles of 
the 5'-matrix. Let us not jump to such conclusions too hastily. Having analyzed how the 
cross section is built up out of individual contributions from the resonances, in the subse- 
quent sections we will see that the cr(n — )■ n') curves are the results of rather complicated 
interplay of all the poles and therefore a Breit-Wigner parametrization of them would be 
naive and misleading. 

6.3 Mittag-Leffler expansion 

The Mittag-Leffler theorem [61] offers a way to decompose the S'-matrix in a sum of terms 
representing its poles and the background. The theorem itself is rather general. In plain 
words, it says that a meromorphic function can be expanded in a sum of an entire function 
and a series of its principal parts at all the poles. What we actually need is more simple. 
Leaving the mathematical subtleties, we want to expand the S'-matrix in a sum of few 
singular terms (poles) and something that takes into account everything else (background 
term). This can be done using the Cauchy integral formula. Indeed, consider a contour 
shown in Fig. [HI which encloses poles of the S-matrix. Then the Cauchy formula reads 




'n 



(69) 



N 



Res[S{Ej)] 



C-E 



dC = 2mS{E) + 2m ^ 



(70) 
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where E is a point inside the contour and Ej are the resonance energies. Therefore in the 
desired expansion, 



N 



SiE) 



^ Res[S{E,)] ^ 



E-E, 



27ii J (-E ^ 



(71) 



the selected poles are represented explicitly while the rest of them as well as the 
background (the Mittag-Leffler's entire function) are taken into account via the contour 
integral. In order to understand what role a particular resonance plays in the scattering 
process, we will simply omit the corresponding term from the sum (ITTll and analyze the 
consequences. In a sense, the expansion (ITTll is similar to that of Ref.jH] where the 
spectral-density function m"'"(A) was presented as a sum over the resonances and the 
background term. It should be noted, however, that Eq. flTTl) is not an approximation 
but an exact expansion when all N poles enclosed by the contour are taken into account. 
To calculate the S'-matrix via Eq. flTTl) . we need to know its residues at the poles and 
to be able to calculate the integral for a given energy E. The residues can be found by 
numerical differentiation of the Jost matrix determinant. Indeed, according to Eq. 



S = f 



(out) 



/(in) 
J 22 

_ f (in) 
J2I 



/(in) 
'J12 



ill 



(in) I det/^i'^) 



(72) 



where the poles are caused by simple zeros of det/('°)(E). Therefore 



Res [S{E)] = /(°"*^(E) 



f22 



d 
dE 



det/('")(E) 



-1 



(73) 



with the derivative numerically obtained within the central difference approximation 



d 

dE 



det/('")(E) 



det /('°)(E + e) - det /(*°)(^ - e) 
Ye 



(74) 



In the calculations, we used e = 10^^ which gave the accuracy of at least 5 digits. This was 
sufficient for the purpose of visual comparing the curves when analyzing the contribution 
to the total cross section from different S'-matrix poles. 

Thus calculated residues of the S'-matrix are given in Table [2] and plotted in Fig. |9l It is 
interesting to notice that they follow regular pattern in the form of anti-clockwise spiral. 
The reason for growth of the spiral radius is that according to Eq. fl50|) 



Res [Snn'm 



r r 
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and the widths become larger every time we go to the next resonance. The reason for 
the anti-clockwise turning of the corresponding phase factor, however, remains unclear. 
As is seen from Fig. [TOl the phases of the residues show the tendency to converge to a 
constant angle, which means that for the far away resonances the points on the spiral 
become closer and closer to each other. 

The contour integral of Eq. flTTl) . enclosing nine [N = 9) resonance points (see Fig. [H]), 
was calculated using 175 points of the Gaussian quadrature formula on each of the three 
line segments that form the contour. With this number of quadrature points we achieved 
four digit accuracy which was more than sufficient for our purpose of visual comparing 
the cross section curves. The vertices of the triangle contour (0.5, 0.5), (0.5,-30.0), and 
(20.1, 0.5) were chosen in such a way that both branching points are outside of it and 
with 9 > 0.257r the contour does not cross the unitary cuts. 

A question may arise: Why did we not include the bound state poles in the contour? There 
are two reasons for this. Firstly, the present work was motivated (as we explained in the 
Introduction) by the controversy in associating the irregularities of the chemical reaction 
cross sections with intermediate resonant states. Therefore, our goal was to develop a 
technique for examining the contributions of the resonance poles, not the bound states. 
Secondly, a contour enclosing the bound states, would require calculating the S-matrix 
via Eq. (H6|l at the points around them, i.e. with Imfc„ > 0. At such points, however, 
the limit (136|) does not exist. In order to reach these points, we would need a negative 
rotation angle with 6 < — 7r/2 which is not possible. There is a way to circumvent this 
difficulty by expanding the solutions j^(™/out) series near E = (see Ref. |1S]), but this 
would require a separate publication. It should be emphasized that the inclusion of only 
nine poles in the contour is not an approximation. The Cauchy theorem f l70|) is exact, 
provided that all enclosed poles are taken into account. All the other poles are taken care 
of by the contour integral. 

6.4 Analysis of the cross sections 

Before starting the analysis, we tested the numerical procedure by comparing the cross 
sections directly obtained from Eqs. ( 13T|) at real energies with ^ = and the corresponding 
cross sections calculated from the expansion fl7ip . The same a{n — )■ n') curves, namely, 
shown in Figs. El El and [71 were obtained in both cases (within four digit accuracy, which 
is visually indistinguishable). Having thus established that the sum f l7T|) reproduces the 
cross sections correctly, we started to investigate the importance of its terms by excluding 
them from the sum and calculating the resulting cross sections. 

First of all, we omitted all the nine poles. To our surprise, the remaining integral part of 
the expansion f lTT]) did not turn out to be a smooth function. Contrary to our expectation. 
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at low energies it gives almost the same structure of peaks and deeps as they are in the 
exact cross sections (see Figs. [HI [T21 and US]). This means that these peaks are of a non- 
resonance origin. They are either the threshold irregularities (known in scattering theory 
as threshold cusps [l5l [62]), or the influence of some other S-matrix poles (especially the 
four bound states that are close to the point E = 0). 

In an attempt to understand what role each individual resonance plays in the scattering 
process, we excluded them one at a time from the sum flTTl) . This procedure can give us 
some indication and hints. In Figs. [TH [151 and [T6] the cross sections a{l — )■ 1), a{l — )■ 2), 
and (t(2 — )■ 2) thus obtained are shown for eight different choices of the resonance to be 
excluded. 

These figures do not show the effect of the exclusion of the first resonance because this 
effect is very simple. When the first resonance is omitted, the narrow deep in cr(l — )■ 2) 
and peak in a{2 2) aX E = 4.768 disappear while the other parts of the cross section 
curves remain unchanged. As far as the elastic cross section a{l — )■ 1) is concerned, the 
resonant cusp ai E = 4.768 (see insert on Fig. [5]) is so weak that it is hardly visible 
and therefore the exclusion of the first resonance pole from the sum f[7T]) is practically 
unnoticeable. The reason for that can be found in Table [5J The Su element of the S- 
matrix has extremely small residue at this pole. As a result the corresponding term in 
the sum (!7T|) for real E is always negligible. Being used to the Breit-Wigner picture, we 
expected to have a sharp energy dependence at around ~ 4.8 in all channels. To our 
surprise, this turned out to be not the case because of smallness of Res(S'ii). 
Apart from the first resonance {E = 4.768), the only place where we can see a clear 
manifestation of a resonance through a distinct peak is the top left corner of Fig. \T5[ 
where one of the peaks in the inelastic cross section a{l — )■ 2) disappears together with 
the second resonance pole while the rest of the curve retains its shape. As far as the 
elastic scattering in the channels 1 and 2 is concerned (see Figs. UM and [T6|) . the only 
information we can deduce from the curves is that the second and the third resonances 
are responsible for the bumps between E = 5 and E = 10, the far resonances do scaling 
of the curves, and no resonances manifest themselves through distinct peaks (except for 
the first one, of course). This is from a general physical view understandable in the sense 
that the fourth and higher resonances are so wide that they, when excited, affect a rather 
wide energy region. 

In general, the cross section analysis based on the expansion (I7T]) . is not straight forward. 
This is not surprising because the contributions of the terms of this expansion into the 
cross section are not linear. Indeed, on the right hand side of Eq. f[69l) in addition to 
squares of these terms there are many interference products whose contributions depend 
on relative phases of the terms. Therefore the full cross section is a result of complicated 
interplay of all resonances and the background term. The analysis of the sum flTTl) itself. 
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i.e. examining the contributions of the resonances into the S'-matrix, could shed some 
hght on the role of individual resonances in the scattering picture. The most convenient 
way of doing this is the analysis of the Argand plots which is done next. 

6.5 Argand plots 

The Argand plots are widely used to depict complex valued functions depending on a real 
parameter. For each value of this parameter, the function value corresponds to a point in 
the complex plane. If the function is continuous, the set of such points forms a curve. 
In the scattering theory, the term Argand plot has a specific meaning, namely, the curve 
in the complex plane, along which the S'-matrix element moves when the (real) collision 
energy increases. Traditionally the Argand plots were mainly used in particle physics. 
There are however many examples of their application in atomic and molecular physics 
as well [631 EH- 

It can be shown (see, for example, Refs. [63] and [65]) that in the absence of resonances the 
phases (arguments) of elements of the S'-matrix monotonously decrease with the increasing 
energy of collision. Therefore these matrix elements as functions of real energy move along 
circular trajectories in the clockwise direction on the Argand plot. Contrary to that, a 
resonance causes the scattering phase shift to increase by vr, which means an increase of 
the S'-matrix phase by 27r (i.e. full anti-clockwise circle). This is, of course, an idealized 
picture for an isolated resonance with weak influence of the background scattering. 
In a vicinity of a resonance, the S'-matrix can always be split in two terms: The background 
and resonance terms, with decreasing and increasing phases, respectively. Their interplay 
may cause the curvature of the Argand trajectory to change. If a resonance is strong, 
the curvature changes so much as making the point on the curve to move along an arc 
in the anti-clockwise direction. With a weak resonance, the curve changes its curvature 
only slightly. 

In Figs. [T71 UHl and dn] the Argand plots for the S-matrix generated by the potential f lBS]) 
are shown in the energy interval 0.5 < E < 20. The dots on the curves mark the integer 
values of the energy, namely, £'=1,2,3,.... As is seen, all three matrix elements have 
three anti-clockwise arcs in the intervals 4.765 < E < 4.775, 6 < < 9, and 9 < E < 20. 
This means that the corresponding peaks aX E = 4.768 and bumps seen between ~ 6 
and ~ 20 in Figs. El El and El are of the resonance origin, despite the fact that not all 
of them are high or sharp. The only question remains: Which poles are responsible for 
these resonances? 

To answer this question, we make use of the expansion fl7T]) . There is no problem of 
interference here because, in contrast to the cross section, the contribution of each pole 
to the S-matrix is linear. Therefore, excluding them one at a time, we can find out which 
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poles the two anti-clockwise arcs of the Argand plots are associated with. 
First of all, as was mentioned before, the first resonance pole does not contribute to 
anything when E moves away from the point E = 4.768. This is because the denominator 
of the corresponding term in the sum flTTl) becomes much lager than the residue. The 
omission of the first resonance term results therefore in the disappearance of the first anti- 
clockwise circles (4.765 < E < 4.775) in Figs. [T71 HSl and [121 The rest of these figures 
remain unaffected. Contrary to that, the second pole makes a noticeable contribution 
everywhere. In Fig. [20] the Argand plots for the three elements of the S'-matrix are shown 
when this pole is excluded from the sum (I7T!) . It is seen that in all the channels there is 
no second anti-clockwise arc between £^ ~ 6 and ~ 9 as compared with Figs. [T71 [TH 
and [ini We therefore can conclude that the anti-clockwise arc in the interval 6 < -E < 9 
is associated with the second pole and the bumps in this energy interval on the cross 
sections shown in Figs. [5l [6l and [TJ are the manifestation of the second resonance pole, 
namely at E = 7.241200 - (V2)1.511912. 

It is seen that with the exclusion of the second resonance pole, we still have anti-clockwise 
arcs for Sii{E), starting from ~ 9, and for 5*21 and 5*22 starting from ~ 6 
and extending to E ~ 20. This means that one or several other resonance poles are 
responsible for their formation. 

In order to examine this, we restore the second pole and exclude the third resonance 
pole from the expansion (jTTl) . The resulting Argand plots are shown in Fig. [21] As was 
expected, the second anti-clockwise arc (6 < < 9) came back together with the second 
pole. The third arc {E > 9) however did not disappear together with the third pole. 
Therefore the third pole is not fully responsible for this arc and thus for the wide bump 
on the cross sections beyond E ^ 9. 

This bump, which is seen in all three Figs. [5] [6] and [7] is a collective effect of several 
resonances, namely, the third, fourth, and so on. With the presence of the third pole, 
the ant i- clockwise trajectories for > 9 on the Argand plots of 5*21 and 5*22 form 
closed loops. When this pole is excluded, the curvature of these loops is reduced to the 
extent that the closed loops become just open arcs. Since they do not disappear entirely, 
we conclude that the fourth and the subsequent poles also contribute to their formation, 
although not so much as the third pole does. 

In order to convince ourselves that the two bumps on the cross sections between ~ 6 
and ~ 20 are of pure resonant nature, we excluded all nine poles from the sum f l7ip . i.e. 
left only the background integral term. The Argand plots thus obtained are given in Fig. 
[22] It is seen that they show no anti-clockwise motion. Only close to = 20 the Argand 
trajectories become flat and perhaps change sign of the curvature. This means that the 
integral term practically does not contribute to formation of the bumps. It takes into 
account all the other poles that are outside the integration contour. The Argand plots 
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given in Fig. |221 show that these outside poles contribute very hghtly at higher energies. 
The flattening of the curves is the indication of this. 

Finally, it is worthwhile to say a few words about the requirement of unitarity which is 
violated when one analyzes the Argand plots using a truncated Mittag-Leffler expansion. 
Since the physical S-matrix is unitary, its matrix elements cannot be outside the unitary 
circle (circle of unit radius) on the Argand plots. This is indeed the case with the curves 
given in Figs. [T71 [TSl and [121 However when we omit some terms in the expansion flTTl) . 
the unitarity is violated. As a result the Argand trajectories are not bound to the unitary 
circle anymore (see Figs. [20] and [21]). This should not be an obstacle. Actually, we are not 
interested in absolute position of the curve in the complex plane. What we look for are 
the segments of anti-clockwise motion. In fact, when in our example all nine resonance 
poles are excluded the unitarity is violated so much that the Argand trajectories shift very 
far from the origin of the complex plane. This is why in Fig. [22] we do not indicate the 
coordinate system (which in this case does not give any meaningful information anyway). 

7 Conclusion 

The method presented in this paper, for solving a two-channel problem is based on first 
order differential equations that are equivalent to the corresponding Schrodinger equation 
but are more convenient for dealing with the resonant phenomena. The main advantage 
of this method is that it enables us to directly calculate the Jost matrix for practically 
any complex value of the energy. The spectral points (bound and resonant states) can 
therefore be located in a rigorous way as zeros of the Jost matrix determinant. 
When calculating the Jost matrix, we solve differential equations and thus at the same 
time obtain the wave function with the correct asymptotic behaviour that is embedded in 
the solution from the outset in an analytical form. The method therefore gives a complete 
solution of the problem in one run. Another advantage of the method used in this paper, 
is that it offers very accurate way of calculating not only total widths of resonances but 
their partial widths as well. 

Although in this paper we deal with a two-channel problem, all the formulae can be easily 
generalized to any number of channels. For this, we only need to change the matrix 
dimensions from 2 to an appropriate number N. 

Charged particles can also be considered within this method with minor modifications. 
For this, we need to use the Coulomb functions Fi{ri, kr) and Ge{ri, kr) instead of the 
Rikkati-Bessel and Riccati-Neumann functions ji{kr) and ne{kr), and their combinations 
Fe{7],kr) ^iGi{ri^kr) instead of the Riccati-Hankel functions h^^\kr) when construct- 
ing the wave functions and deriving the equations. For — )■ such equations will be 
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automatically transformed into the equations given here. 

The possibility of calculating the Jost matrix and the S'-matrix for any complex energy in 
the resonance domain, enables us not only to locate the S'-matrix resonance poles but to 
calculate its residues at these poles. This, in turn, makes it possible to obtain the Mittag- 
LefHer type expansion of the S-matrix as a sum of the singular terms (representing the 
resonances) and the background term (contour integral). Such an expansion is very useful 
in analyzing the contribution of individual resonances into the scattering. 
The analysis of the cross section turned out to be less informative than expected. The 
reason for this is that the cross section being a quadratic function of the S'-matrix, involves 
all possible interference terms which smear the contributions from individual resonances 
and make it difficult to clearly separate them. However, using the Argand plot technique, 
one can, as demonstrated above, elucidate the origin of the structures in a computed cross 
section as being the effects of a single isolated resonance, a collective contribution of a 
set of resonances, or a non-resonant background scattering (which in fact is a collective 
contribution from all far poles of the S-matrix). 

The two-channel model by Noro and Taylor, considered in this paper as an example, 
exhibits these difficulties in an extreme form. This is because its spectrum involves a 
string of almost completely overlapping resonances. Even in such an extreme case, with 
the proposed approach, we managed to assign the irregularities of the cross section to the 
resonances and the background term. 

Although this model might look artificial, the spectrum of resonances generated by the 
Noro and Taylor potential, is similar to the spectra of some real atomic systems. For 
example, the existence of a string of overlapping resonances that starts with a narrow 
one, is a feature that characterizes the charge-transfer reaction N^~^ + if — >■ iV^"*" -\- 
[T5l [m [T7]. The analysis of the Noro and Taylor model, given in the present paper, could 
therefore shed some light on the problems associated with realistic systems. Furthermore, 
with small numerical modifications, the methods presented here can be directly applied 
to this charge-transfer reaction. The contribution from each resonant eigenstate can 
then be quantified by using the corresponding S'-matrix residue and the Argand plots as 
demonstrated in this paper. 

Such an analysis, may eventually, contribute to the understanding of the FHD resonances 
[3 Uni HI] discussed in the Introduction. Finally, the complex interference between dif- 
ferent resonant states in building the observable quantity - the cross section - tells us 
that it is absolutely necessary to have theoretical tool, like the present, in order to an- 
alyze the high-resolution, low-energy atom-atom, atom-ion and atom-molecule collision 
experiments that eventually may be performed [29] . 
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Table 1: Spectral points generated by the potential (168|) and shown in Fig. HI 
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Table 2: Residues of the S*- matrix elements at the 9 resonance poles enclosed in the 
contour shown in Fig. [8l 
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Figure 1: Fragments of the physical (Si) and unphysical {S2) sheets of the complex-energy 
Riemann surface around a branching point {B) that corresponds to a threshold energy [E2 or 
£^1). Transition from to S2 and back is possible through the unitary cut running from B 
to infinity. 



ImE 




Figure 2: Typical distribution of the bound states and resonances i.e. spectral points (open 
circles) in the complex energy-plane. It is assumed that Ei < E2 and the energy is measured 
relative to the lowest threshold Ei. The unitary cuts going from the branching points to infinity 
are also shown. Because of the complex rotation, these cuts are turned into the unphysical 
sheet by the angle 29. 
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Figure 3: The integration path for solving Eqs. (I66l) on the segment ob, and the Jost 
function equations f l3T]) on 6/2. 
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Figure 4: Spectral points generated by the potential (1681) and given in Table [H 
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Figure 5: Energy dependence of the elastic scattering cross section in channel 1 for the 
potential (1551) . Few of the S-matrix poles (see Tableland Fig.SD are shown in the lower part 
of the Figure. 
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Figure 6: Cross section energy dependence of the inelastic transition (1 — 2) for the potential 
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Figure 7: Energy dependence of the elastic scattering cross section in channel 2 for the 
potential ( l68i l. 
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Figure 8: The integration contour for Eq. ( iTOl) . The vertices of the triangle are (0.5,0.5), 
(0.5,-30), and (20.1,0.5). 
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Figure 9: Residues of the ^-matrix elements at the 9 resonance poles enclosed in the contour 
shown in Fig. |H1 The coordinates of the points are given in Table [21 
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Figure 10: Phases of the residues of the 5'-matrix elements at the first 14 resonance poles 
given in Table [H 
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5 10 15 

Figure 11: Comparison of the exact cross section for the elastic scattering in the channel 1 
(thin curve) with the corresponding cross section obtained from the expansion (ITTl) where all 
the pole terms are omitted (thick curve). 




Figure 12: Comparison of the exact cross section for the inelastic transition (1 — )■ 2) (thin 
curve) with the corresponding cross section obtained from the expansion ( ITTl) where all the 
pole terms are omitted. In order to fit into the picture the thick curve is scaled down by the 
factor of 1/100. 
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Figure 13: Comparison of the exact cross section for the elastic scattering in the channel 2 
(thin curve) with the corresponding cross section obtained from the expansion (171 ft where all 
the pole terms are omitted. In order to fit into the picture the thick curve is scaled down by 
the factor of 1/10. 
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Figure 14: Exact cross section for the elastic scattering in channel 1 (thin curve 
corresponding cross section obtained from the expansion (TTTD where one of the 
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Figure 15: Exact cross section for the inelastic transition 1—7-2 (thin curve) and the cor- 
responding cross section obtained from the expansion (ITTI) where one of the nine pole terms 
is omitted (thick curve). As indicated on thg graphs, some of the thick curves are scaled 
down by the factor of 0.1 in order to fit into the picture. Vertical dashed lines indicate the 
positions Ej. of the excluded resonances and thick horizontal bars on the energy axis cover 
the corresponding intervals Ej- ± T/2. For the last four resonances, such intervals exceed the 
energy segment shown on the figure. 
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Figure 16: Exact cross section for the elastic scattering in channel 2 (thin curve) and the 
corresponding cross section obtained from the expansion (17 ID where one of the nine pole 
terms is omitted (thick curve). As indicated^gpn the graphs, some of the thick curves are 
scaled down by the factor of 0.1 in order to fit into the picture. Vertical dashed lines indicate 
the positions of the excluded resonances and thick horizontal bars on the energy axis cover 
the corresponding intervals Ej- ± T/2. For the last four resonances, such intervals exceed the 
energy segment shown on the figure. 
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Figure 17: Argand plot of Sii{E) in the energy interval from Ei = 0.5 to E2 = 20. The dots 
on the curve indicate the points corresponding to E = 1, 2, 3, 
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Figure 18: Argand plot of S2i{E) in the energy interval from Ei = 0.5 to E2 = 20. The dots 
on the curve indicate the points corresponding to £^ = 1, 2, 3, 
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Figure 20: Argand plot of the S'-matrix with the second pole excluded. The dots mark the 
points corresponding to E = 6 and E = 9. 




Figure 21: Argand plot of the ^-matrix with the third pole excluded. The dots nnark the 
points corresponding to £^ = 6 and E — 9. 
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Figure 22: Argand plot of the ^'-matrix with all nine resonance poles excluded. 
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